*****************************************************************************************************
* DIVERSIONARY DESECRATION? Regime Instability and Societal Violence against Minority Sacred Spaces *
* Ariel Zellman and Andrea Malji; Politics, Religion & Ideology (2023) *
************************************************************************

******************************
* Site Targeting Main Models *
******************************

* MODEL 1 *

mlogit wzsocdis_sitetargeting_alt mmxx lxx wzmin2majx v2x_polyarchy cpi_pchg unempl_youthm_ilo max_couprisk zpop ezmajpctx minpctx fraction_rasm polarize_rasm shared_religion i.year, vce(cluster group_cow) nolog

* Model 1 Marginal Effects *

margins, dydx(mmxx wzmin2majx cpi_pchg unempl_youthm_ilo max_couprisk v2x_polyarchy lxx) predict(outcome(1))

margins, dydx(mmxx wzmin2majx cpi_pchg unempl_youthm_ilo max_couprisk v2x_polyarchy lxx) predict(outcome(2))

margins, dydx(mmxx wzmin2majx cpi_pchg unempl_youthm_ilo max_couprisk v2x_polyarchy lxx) predict(outcome(3))

* Model 1 Marginal Probabilities *

margins, at(mmxx=(0 2 24)) predict(outcome(1))

margins, at(mmxx=(0 2 24)) predict(outcome(2))

margins, at(mmxx=(0 2 24)) predict(outcome(3))

margins, at(wzmin2majx=(0 3 6)) predict(outcome(1))

margins, at(wzmin2majx=(0 3 6)) predict(outcome(2))

margins, at(wzmin2majx=(0 3 6)) predict(outcome(3))

margins, at(cpi_pchg=(0.060 4.827 43.073)) predict(outcome(1))

margins, at(cpi_pchg=(0.060 4.827 43.073)) predict(outcome(2))

margins, at(cpi_pchg=(0.060 4.827 43.073)) predict(outcome(3))

margins, at(unempl_youthm_ilo=(2.425 15.803 39.835)) predict(outcome(1))

margins, at(unempl_youthm_ilo=(2.425 15.803 39.835)) predict(outcome(2))

margins, at(unempl_youthm_ilo=(2.425 15.803 39.835)) predict(outcome(3))

margins, at(max_couprisk=(0.00003 0.0017 0.0119)) predict(outcome(1))

margins, at(max_couprisk=(0.00003 0.0017 0.0119)) predict(outcome(2))

margins, at(max_couprisk=(0.00003 0.0017 0.0119)) predict(outcome(3))

margins, at(v2x_polyarchy=(0.095 0.532 0.892)) predict(outcome(1))

margins, at(v2x_polyarchy=(0.095 0.532 0.892)) predict(outcome(2))

margins, at(v2x_polyarchy=(0.095 0.532 0.892)) predict(outcome(3))

margins, at(lxx=(3 7 29)) predict(outcome(1))

margins, at(lxx=(3 7 29)) predict(outcome(2))

margins, at(lxx=(3 7 29)) predict(outcome(3))

* MODEL 2 - Violence v Sacred Sites *

mlogit wzsocdis25x mmxx lxx wzmin2majx v2x_polyarchy cpi_pchg unempl_youthm_ilo max_couprisk zpop ezmajpctx minpctx fraction_rasm polarize_rasm shared_religion i.year, vce(cluster group_cow)

* Model 2 Marginal Effects, Violence *

margins, dydx(mmxx wzmin2majx cpi_pchg unempl_youthm_ilo max_couprisk v2x_polyarchy lxx) predict(outcome(3))

* Model 2 Marginal Probabilities, Violence *

margins, at(mmxx=(0 2 24)) predict(outcome(3))

margins, at(wzmin2majx=(0 3 6)) predict(outcome(3))

margins, at(cpi_pchg=(0.060 4.827 43.073)) predict(outcome(3))

margins, at(unempl_youthm_ilo=(2.425 15.803 39.835)) predict(outcome(3))

margins, at(max_couprisk=(0.00003 0.0017 0.0119)) predict(outcome(3))

margins, at(v2x_polyarchy=(0.095 0.532 0.892)) predict(outcome(3))

margins, at(lxx=(3 7 29)) predict(outcome(3))

* MODEL 3 - Violence v Minorities *

mlogit wzsocdis23x mmxx lxx wzmin2majx v2x_polyarchy cpi_pchg unempl_youthm_ilo max_couprisk zpop ezmajpctx minpctx fraction_rasm polarize_rasm shared_religion i.year, vce(cluster group_cow)

* Model 3 Marginal Effects, Violence *

margins, dydx(mmxx wzmin2majx cpi_pchg unempl_youthm_ilo max_couprisk v2x_polyarchy lxx) predict(outcome(3))

* Model 3 Marginal Probabilities, Violence *

margins, at(mmxx=(0 2 24)) predict(outcome(3))

margins, at(wzmin2majx=(0 3 6)) predict(outcome(3))

margins, at(cpi_pchg=(0.060 4.827 43.073)) predict(outcome(3))

margins, at(unempl_youthm_ilo=(2.425 15.803 39.835)) predict(outcome(3))

margins, at(max_couprisk=(0.00003 0.0017 0.0119)) predict(outcome(3))

margins, at(v2x_polyarchy=(0.095 0.532 0.892)) predict(outcome(3))

margins, at(lxx=(3 7 29)) predict(outcome(3))

******************************************************
* Site Targeting Paper, Robustness Checks on Model 1 *
******************************************************

* Substitute Polity for VDem

mlogit wzsocdis_sitetargeting_alt mmxx lxx wzmin2majx zpolity2 cpi_pchg unempl_youthm_ilo max_couprisk zpop ezmajpctx minpctx fraction_rasm polarize_rasm shared_religion i.year, vce(cluster group_cow) nolog

* Substitute Polity and Polity-squared for VDem

mlogit wzsocdis_sitetargeting_alt mmxx lxx wzmin2majx zpolity2 zpolity2_sq cpi_pchg unempl_youthm_ilo max_couprisk zpop ezmajpctx minpctx fraction_rasm polarize_rasm shared_religion i.year, vce(cluster group_cow) nolog

* Sub CoupRisk for CNTS weighted internal conflict measure

mlogit wzsocdis_sitetargeting_alt mmxx lxx wzmin2majx v2x_polyarchy cpi_pchg unempl_youthm_ilo cnts_conflict_wt zpop ezmajpctx minpctx fraction_rasm polarize_rasm shared_religion i.year, vce(cluster group_cow) nolog

* Sub CoupRisk for CNTS elite unrest measure

mlogit wzsocdis_sitetargeting_alt mmxx lxx wzmin2majx v2x_polyarchy cpi_pchg unempl_youthm_ilo unrest_elite zpop ezmajpctx minpctx fraction_rasm polarize_rasm shared_religion i.year, vce(cluster group_cow) nolog

* Control for Violence v RelMin

mlogit wzsocdis_sitetargeting_alt wzsocdis23x mmxx lxx wzmin2majx v2x_polyarchy cpi_pchg unempl_youthm_ilo max_couprisk zpop ezmajpctx minpctx fraction_rasm polarize_rasm shared_religion i.year, vce(cluster group_cow) nolog

* Alternative Model 1 using 4-level Site Targeting Variable, Major Vandalism and Minor Violence as Distinct Levels

mlogit wzsocdis_sitetargeting4 mmxx lxx wzmin2majx v2x_polyarchy cpi_pchg unempl_youthm_ilo max_couprisk zpop ezmajpctx minpctx fraction_rasm polarize_rasm shared_religion i.year, vce(cluster group_cow) nolog

margins, dydx(mmxx lxx wzmin2majx cpi_pchg unempl_youthm_ilo max_couprisk v2x_polyarchy) predict(outcome(1))

margins, dydx(mmxx lxx wzmin2majx cpi_pchg unempl_youthm_ilo max_couprisk v2x_polyarchy) predict(outcome(2))

margins, dydx(mmxx lxx wzmin2majx cpi_pchg unempl_youthm_ilo max_couprisk v2x_polyarchy) predict(outcome(3))

margins, dydx(mmxx lxx wzmin2majx cpi_pchg unempl_youthm_ilo max_couprisk v2x_polyarchy) predict(outcome(4))
